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ABSTRACT 



We present the scaling relation between Sunyaev-Zeldovich (SZ) signal and stellar mass for almost 260,000 locally brightest galaxies (LBGs) 
selected from the Sloan Digital Sky Survey (SDSS). These are predominantly the central galaxies of their dark matter halos. We calibrate the 
stellar-to-halo mass conversion using realistic mock catalogues based on the Millennium Simulation. Applying a multi-frequency matched filter to 
the Planck data for each LBG, and averaging the results in bins of stellar mass, we measure the mean SZ signal down to ~ 2 x 10^^ M©, with a 
clear indication of signal at even lower stellar mass. We derive the scaling relation between SZ signal and halo mass by assigning halo properties 
from our mock catalogues to the real LBGs and simulating the Planck observation process. This relation shows no evidence for deviation from 
\ a power law over a halo mass range extending from rich clusters down to M500 ~ 2 x 10^^ M©, and there is a clear indication of signal down to 
. ^ ' ^500 ~ 4 X 10^^ Mq. Planck's SZ detections in such low-mass halos imply that about a quarter of all baryons have now been seen in the form of hot 
halo gas, and that this gas must be less concentrated than the dark matter in such halos in order to remain consistent with X-ray observations. At 
the high-mass end, the measured SZ signal is 20 % lower than found from observations of X-ray clusters, a diff'erence consistent with Malmquist 
^ bias eff'ects in the X-ray sample. 

Key words, cosmology: observations — cosmic microwave background — large-scale structure of the Universe — galaxies: clusters: general 



1. Introduction 

Galaxy evolution is currently understood to reflect a thermal cy- 
cle operating between baryonic components confined in dark 
matter halos. Gas cools radiatively during the hierarchical build- 
up of the halo population and condenses to form galaxies in halo 
cores. Left unchecked, cooling results in more massive galaxies 
than observed (Balogh et al. 2001; Lin & Mohr 2004; Tomatore 



* Corresponding author: J. A. Rubino-Martm, jalberto@iac . es 



et al. 2003), and one must invoke an additional source of non- 
gravitational heating to prevent a "cooling crisis" (White & Rees 
1978; Cole 1991; White & Frenk 1991; Blanchard et al. 1992). 
Feedback from star formation and supernovae appears insufift- 
cient to halt cooling in massive halos (Borgani et al. 2004), so 
some modelers have invoked additional heating by active galac- 
tic nuclei (AGN, Churazov et al. 2002; Springel et al. 2005a; 
McNamara & Nulsen 2007). Such models show substantially 
improved agreement with the luminosity-temperature relation 
of X-ray clusters (Valageas & Silk 1999; Bower et al. 2001; 
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Cavaliere et al. 2002) and the luminosity function of galaxies 
(Croton et al. 2006; Bower et al. 2006; Somerville et al. 2008). 
The energetics of AGN feedback imply that it should have es- 
pecially strong effects on low-mass clusters, heating gas in the 
central regions and pushing it to larger radii, thereby reducing 
both gas fractions and X-ray luminosities (Puchwein et al. 2008; 
McCarthy et al. 2010). 

Relationships between the gas, stellar, and dark matter prop- 
erties of halos are important to our understanding of galaxy for- 
mation. Measurements of these relationships over a wide range 
of halo mass, from rich clusters down to individual galaxies, are 
therefore a primary objective of a number of current observa- 
tional campaigns. Recent studies have probed the relationship 
between the mass of a halo and the stellar mass of its central 
galaxy (the SHM relation) using "abundance matching" tech- 
niques, the dynamics of satellite galaxy populations, and gravita- 
tional lensing (Guo et al. 2010; Moster et al. 2010; Mandelbaum 
et al. 2006; Leauthaud et al. 201 1). 

Corresponding constraints on the gas content of halos over 
a similar mass range are not yet available. Although there are 
many detailed X-ray studies of the intracluster medium, these 
mostly concern massive clusters; lower mass groups are faint 
and so are difficult to study individually. The Sunyaev-Zeldovich 
(SZ) effect (Sunyaev & Zeldovich 1972; Birkinshaw 1999) of- 
fers a fresh means to address this problem. Large-area SZ sur- 
veys are just beginning to be amassed by ground-based instru- 
ments such as the Atacama Cosmology Telescope (ACT, Swetz 
et al. 2008; Marriage et al. 2010; Sehgal et al. 2010; Hand et al. 
2011), the South Pole Telescope (SPT, Carlstrom et al. 2009; 
Staniszewski et al. 2009; Vanderlinde et al. 2010; Williamson 
et al. 2011) and APEX-SZ (Dobbs et al. 2006), as well as by 
the Planck^ satellite mission, (Planck Collaboration VIII 2011; 
Planck Collaboration IX 2011; Planck Collaboration X 2011; 
Planck Collaboration XI 201 1 ; Planck Collaboration XII 201 1). 

High S/N observations of individual objects are not currently 
possible over the full mass range from galaxy clusters down to 
individual bright galaxies. The SHM relation can only be esti- 
mated for lower mass objects through statistical methods applied 
to large catalogues. In this context, the SZ effect presents excit- 
ing new opportunities. First steps in this direction were taken by 
Planck Collaboration XII (2011) and Hand et al. (2011), with 
more recent work by Draper et al. (2012) and Sehgal et al. 
(2012). In our first study (Planck Collaboration XII 2011), we 
binned large numbers of maxBCG (Koester et al. 2007) clusters 
by richness to measure the relation between mean SZ signal and 
richness. In a similar manner. Hand et al. (2011) binned ACT 
measurements of luminous red galaxies to determine the mean 
relation between SZ signal and LRG luminosity. 

Here, we extend our previous work with Planck multi- 
frequency observations of a large sample of locally brightest 
galaxies (LBGs). These were selected from the Sloan Digital 
Sky Survey (SDSS) using criteria designed to maximize the frac- 
tion of objects that are the central galaxies of their dark matter 
halos. We stack the Planck data in order to estimate the mean 
SZ signal for LBGs in a series of stellar mass bins. We then 
use mock galaxy catalogues based on the Millennium Simulation 
and tuned to fit the observed abundance and clustering of SDSS 



^ Planck (http://www.esa.int/Planck) is a project of the 
European Space Agency (ESA) with instruments provided by two sci- 
entific consortia funded by ESA member states (in particular the lead 
countries France and Italy), with contributions from NASA (USA) and 
telescope reflectors provided by a collaboration between ESA and a sci- 
entific consortium led and funded by Denmark. 



galaxies to establish the relation between stellar and halo mass. 
Planck is a unique SZ instrument for this purpose because of its 
large frequency coverage and the fact that it observes the entire 
SDSS survey area, allowing study of large samples of galaxy 
systems with extensive multi- wavelength data. 

We unambiguously (> 3cr) detect the SZ signal down to stel- 
lar masses of 2 x 10^^ M©, corresponding to an effective halo 
mass M500 of 2 X 10^^ M© (see Sect. 2) and we find clear indica- 
tions of signal down to 10^^ M© (M500 = 4 x 10^^ M©). Detailed 
simulation both of the galaxy sample and of the Planck mea- 
surement process allows us to correct the effects of halo mis- 
centering and of the scatter in halo mass at fixed stellar mass 
when estimating the SZ signal-halo mass relation. We find that 
the relation is well described by a single power law within its 
statistical uncertainties. At the high end, our results overlap the 
mass range probed by X-ray clusters, where we find a 20% lower 
SZ signal than obtained from fits to X-ray selected cluster sam- 
ples. This difference is consistent with possible Malmquist bias 
effects in the X-ray sample. The gas properties of dark matter 
halos display a remarkable regularity from the poorest groups to 
the richest clusters. 

Throughout this paper, we adopt a fiducial ACDM cosmol- 
ogy consistent with the WMAP7 results (Komatsu et al. 2011). 
In particular, we use = 0.272, Qa = 0.728, = 0.961, 
and CTg = 0.807. We express the Hubble parameter at redshift 
z as H(z) = HoE(z), with Hq = h x lOOkms'^ Mpc"^ and 
h = 0.704. For the redshift range of interest (z ^ 1), we ap- 
proximate E^(z) = ^^m(l + z)^ + ^^A. The virial radius of a halo 
is defined here as R200, the radius enclosing a mean density 200 
times the critical density at that redshift, i.e., 200 x pdz), where 
Pc(z) = 3//^(z)/(8;rG). The virial mass is then defined as 

M200 = 200(47r/3)7?^ooPe, 

which we also refer to as Mh. Similarly, we quote the conven- 
tional masses M500 and radii, R500, when presenting the SZ scal- 
ings. For stellar mass, we use the symbol M*. 

The SZ signal is characterized by Ysoo, the Comptonization 
parameter integrated over a sphere of radius R500, expressed in 
square arcminutes. Specifically, 

Ysoo^((TT/(m,c^)) PdV/Dl(z), 
Jo 

where Da(z) is the angular-diameter distance, crj is the Thomson 
cross-section, c is the speed of light, Mq is the electron rest 
mass, and P = n^kT^ is the pressure, obtained as the product 
of the electron number density and the electron temperature. 
Throughout this paper, we use the quantity 

?500 = YsmE-^l\z){DA{z)l5mM^c)\ 

also expressed in square arcminutes, as the intrinsic SZ signal, 
scaled to redshift z = and to a fixed angular diameter distance. 

The paper is organized as follows. Sect. 2 describes the 
Planck maps used in our analysis, and our reference catalogue 
of locally brightest galaxies, based on SDSS data. Sect. 3 de- 
scribes our methodology. Sects. 4 and 5 gives our main results 
and the tests made to demonstrate their robustness. Sections 6 
and 7 contain discussion and conclusions, respectively. 

2. Data 

2.1. Planck Data Set 

Planck (Tauber et al. 2010; Planck Collaboration I 2011) is the 
third generation space mission to measure the anisotropy of the 
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cosmic microwave background (CMB). It observes the sky in 
nine frequency bands covering 30-857 GHz with high sensitiv- 
ity and angular resolution from 3r to 5'. The Low Frequency 
Instrument (LFI; Mandolesi et al. 2010; BersanelH et al. 2010; 
Mennella et al. 201 1) covers the 30, 44, and 70 GHz bands with 
amplifiers cooled to 20 K. The High Frequency Instrument (HFI; 
Lamarre et al. 2010; Planck HFI Core Team 2011a) covers the 
100, 143, 217, 353, 545, and 857 GHz bands with bolometers 
cooled to 0.1 K. Polarisation is measured in all but the highest 
two bands (Leahy et al. 2010; Rosset et al. 2010). A combina- 
tion of radiative cooling and three mechanical coolers produces 
the temperatures needed for the detectors and optics (Planck 
Collaboration II 2011). Two data processing centres (DPCs) 
check and calibrate the data and make maps of the sky (Planck 
HFI Core Team 201 lb; Zacchei et al. 201 1). Planck'^ sensitivity, 
angular resolution, and frequency coverage make it a powerful 
instrument for Galactic and extragalactic astrophysics as well 
as for cosmology. Early astrophysics results are given in Planck 
Collaboration VIII-XXVI 2011, based on data taken between 
13 August 2009 and 7 June 2010. Intermediate astrophysics re- 
sults are now being presented in a series of papers based on data 
taken between 13 August 2009 and 27 November 2010. 

2.2. A Locally Brightest Galaxy Catalogue 

To select a sample of central galaxies, we first define a par- 
ent population with r < 17.7 (r-band, extinction-corrected, 
Petrosian magnitude) from the spectroscopic galaxy catalogue 
of the New York University Value Added Galaxy Catalogue^. 
This was built by Blanton et al. (2005) based on the seventh data 
release of the Sloan Digital Sky Survey (SDSS/DR7 Abazajian 
et al. 2009). This parent catalogue contains 602,251 galaxies. 
We then define "locally brightest galaxies" to be the set of all 
galaxies with z > 0.03 that are brighter in r than all other sample 
galaxies projected within 1.0 Mpc and with redshift diff'ering by 
less than 1,000km s"^. After this cut 347,486 locally brightest 
galaxies remain. 

The SDSS spectroscopic sample is incomplete because it 
proved impossible to place a fibre on every object satisfying the 
photometric selection criteria, and because some spectra failed 
to give acceptable redshifts. The completeness to our chosen 
magnitude limit varies with position, with a mean of 91 .5 % over 
the survey as a whole. To ensure that galaxies without SDSS 
spectroscopy do not violate our sample selection criteria, we 
have used SDSS photometry to eliminate all objects with a com- 
panion that is close and bright enough that it might violate the 
above criteria. Specifically, we have used the "photometric red- 
shift 2" catalogue (photoz2 Cunha et al. 2009) from the SDSS 
DR7 website to search for additional companions. This cata- 
logue tabulates a redshift probability distribution in bins of width 
Az = 0.0145 for every galaxy down to photometric limits much 
fainter than we require. We then eliminate any candidate with 
a companion in this catalogue of equal or brighter r-magnitude 
and projected within 1.0 Mpc, unless the photometric redshift 
distribution of the "companion" is inconsistent with the spectro- 
scopic redshift of the candidate. (Our definition of "inconsistent" 
is that the total probability for the companion to have a redshift 
equal to or less than that of the candidate is less than 0.1; in 
practice this eliminates "companions" that are too red to be at a 
redshift as low as that of the candidate.) This procedure leaves 
us with a cleaned sample of 259,579 locally brightest galaxies. 
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The NYU-VAGC provides a variety of data for each galaxy. 
In addition to the positions, magnitudes, and redshifts used to 
create our sample, we will make use of rest-frame colours and 
stellar masses. The latter are based on stellar population fits to 
the five-band SDSS photometry and on the measured redshifts, 
assuming a Chabrier (2003) stellar initial mass function (Blanton 
& Roweis 2007). In Fig. 1 we compare the colour and redshift 
distributions of our final sample of locally brightest galaxies to 
those of the parent sample for five disjoint ranges of stellar mass. 
For logiQ M^/Mq > 10.8, the distributions are similar for the 
two populations. At lower stellar mass, locally brightest galaxies 
are a small fraction of the parent sample and are biased to bluer 
colours and to slightly larger redshifts. In our stacking analysis 
below, we obtain significant SZ signals only for galaxies with 
logiQ M:,/Mo > 11.0. Our sample contains 81,392 galaxies sat- 
isfying this bound, the great majority of them on or near the red 
sequence. 

2.2.1 . The reliability of our central galaxy sample and its 
stellar mass-halo mass relation 

We expect the majority of our locally brightest galaxies to be 
the central galaxies of their dark matter halos, just as bright field 
galaxies lie at the centres of their satellite systems and cD galax- 
ies lie near the centres of their clusters and are normally their 
brightest galaxies. For our later analysis, it is important to know 
both the reliability of our galaxy sample, i.e., the fraction of 
galaxies that are indeed the central galaxies of their halos, and 
the relation between the observable stellar masses of the galax- 
ies and the unobservable masses of their halos. In this section we 
investigate both issues using an update of the publicly available^ 
semi-analytic galaxy formation simulation of Guo et al. (2011). 
The update uses the technique of Angulo & White (2010) to 
rescale the Millennium Simulation (Springel et al. 2005b) to the 
WMAP7 cosmology, then readjusts the galaxy formation param- 
eters to produce a z = galaxy population with abundance and 
clustering properties that are almost indistinguishable from those 
of the original model. At the relatively high masses relevant for 
our work, this simulation provides a very close match to the ob- 
served luminosity and stellar mass functions of the SDSS as well 
as to the auto-correlations of SDSS galaxies as a function of stel- 
lar mass (Guo et al. 2012). 

We construct a sample of locally brightest galaxies from this 
simulation using criteria exactly analogous to those used for the 
measured data. We project the galaxy distribution onto one of 
the faces of the simulation cube and assign each galaxy a red- 
shift based on its distance and peculiar velocity in the projection 
direction. A galaxy is considered locally brightest if it has no 
neighbour that is brighter in r within 1 .0 Mpc projected distance 
and 1,000 km s"^ in redshift. We divide galaxies into "centrals", 
defined as those lying at the minimum of the gravitational poten- 
tial of the dark matter friends-of-friends (FoF) group with which 
they are associated, and "satellites", defined as all other galaxies. 

With these definitions we can assess the fraction of our lo- 
cally brightest galaxies that are truly central galaxies. The black 
line in Fig. 2 shows, as a function of stellar mass, the fraction of 
all galaxies in the simulation that are centrals. At stellar masses 
just above 10^^ M© this fraction is about one half, but it increases 
with stellar mass, reaching two thirds by log^QM^^/M© = 11.0 
and 90% by log^QM^^/M© = 11.8. In contrast, the fraction of 
locally brightest galaxies that are centrals is much higher, with 
a minimum of just over 83 % at stellar masses somewhat above 

^ http : //www . mpa- garching . mpg . de/millennium 
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Fig. 1. Distributions in colour (left) and redshift (right) of our locally brightest galaxies and of the SDSS/DR7 population from 
which they were drawn. Black histograms refer to the parent sample and red histograms to the locally brightest galaxies. The panels 
in each set correspond to five disjoint ranges of log^Q M^/Mq, as indicated in the labels. In the left-hand set, additional labels give 
the number of galaxies contributing to the parent (black) and locally brightest (red) histograms. Dashed vertical lines in these same 
panels indicate the colour we use to separate red and blue galaxies in Fig. 3 below. 
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Fig. 2. Fraction of locally brightest galaxies that are the central 
objects in their dark halos, based on the simulations of Guo et al. 
(201 1). The solid line traces the fraction of all simulated galaxies 
that are central galaxies as a function of stellar mass. This frac- 
tion increases with stellar mass, reaching 90 % at the high mass 
end. The dashed line presents the central galaxy fraction for lo- 
cally brightest galaxies selected from the simulation according 
to the criteria applied to the SDSS data. This yields a sample 
that is over 83 % reliable at all stellar masses. 



10^^ Mq. We have checked those locally brightest galaxies that 
are satellites, finding that for log^QM^/M© > 11, about two- 
thirds are brighter than the true central galaxies of their halos. 
The remainder are fainter than their centrals, and are considered 
locally brightest because they are more than 1 Mpc (projected) 
from their centrals (60%) or have redshifts diff'ering by more 
than 1,000km s-i (40%). 



We can assign a halo mass, M200, to every galaxy in our sim- 
ulation. For both satellite galaxies and central galaxies, we take 
M200 to be the current M200 of the FoF dark matter halo with 
which the object is associated, i.e., the mass contained within 
its 7?2oo- Figure 3 shows a scatter plot of M200 against for a 
random subset (one out of every 80) of our sample of simulated 
locally brightest galaxies. We indicate central galaxies with red 
or blue points according to their rest-frame g - r colour (with 
the two distinct regions separated by the vertical dashed lines in 
the left panel of Fig. 1) while satellite galaxies are indicated by 
black points. Clearly, red (passive) and blue (star-forming) cen- 
tral galaxies lie on diff'erent M2oo-^* relations. That for passive 
galaxies is steeper, and is off'set to larger halo mass in the stellar 
mass range where both types of central galaxy are present. 

Satellite galaxies lie in halos in the massive tail of the dis- 
tribution for central galaxies of the same stellar mass. Satellites 
misidentified as centrals in our catalogue are usually outlying 
members projected at relatively large separation (from a few 
hundred kiloparsecs to 2 Mpc). Their presence bias high both the 
mean halo mass (the high black points in Fig. 3) and the spatial 
extent of the stacked SZ signals we measure below. However, 
since two thirds of the satellites that we misidentify as central 
galaxies are in fact brighter than the true central galaxies of their 
halos (i.e., they are not typical satellites), this bias is not extreme. 
In any case, we correct for these eff'ects explicitly in our analysis 
using the simulation. 

The lower of the two continuous curves in Fig. 3 shows the 
median M200 as a function of M^. We will take this as an esti- 
mate of the typical halo mass associated with a central galaxy of 
known M^, and will use it to set the angular size of the matched 
filter for each observed galaxy when stacking SZ signal as a 
function of stellar mass. The upper continuous curve shows the 
mean M200 as a function of . The substantial shift between the 
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Fig. 3. Scatter plot of M200 against for a random subset (one 
out of 80) of our sample of simulated locally brightest galaxies. 
Central galaxies are shown as red or blue points according to 
their g - r colour, using the cuts indicated in Fig. 1. Satellite 
galaxies are shown as black points. The lower and upper curves 
give the median and mean values of halo mass as a function of 
stellar mass. 



two is a measure of the skewness induced by the differing rela- 
tions for passive and star-forming centrals and by the presence 
of the tail of cluster satellite galaxies (see Appendix B). 

3. Analysis 

Our analysis closely follows that presented in Planck 
Collaboration X (2011), Planck Collaboration XI (2011), and 
Planck Collaboration XII (201 1), employing as primary method 
a multi-frequency matched filter (hereafter MMF) optimized in 
both frequency and angular space to extract the thermal SZ sig- 
nal (Herranz et al. 2002; Melin et al. 2006). We find that dust 
emission from our target sources aff'ects the MMF measurements 
noticeably at low stellar mass, and that an eff'ective mitigation is 
to restrict our final measurements to the three lowest HFI fre- 
quencies (100, 143, and 217 GHz). This is detailed in Sect. 5. 
Our primary scientific results are hence all based on this three- 
band MMF. 

For the SZ model template, we employ, as in earlier work 
(Planck Collaboration X 2011; Planck Collaboration XI 2011; 
Planck Collaboration XII 2011), the so-called "universal pres- 
sure profile" (Arnaud et al. 2010) deduced from X-ray observa- 
tions of the REXCESS cluster sample (Bohringer et al. 2007). 
The /?5oo value associated with the halo of each central galaxy 
is obtained as follows. We first use the SHM relation giving the 
median halo M200 as a function of central galaxy stellar mass, as 
presented in Sect. 2.2.1. Then, using an NEW profile (Navarro 
et al. 1997) and the concentration parameter C200 given by Neto 
et al. (2007), we convert M200 to M500 and derive R^m for each 
halo. The angular scale for the filter is finally given by projecting 
7?5oo at the redshift of the target LBG. 

In addition to the MMF, and in order to test the robustness of 
the results, the impact of foreground contamination and possible 



systematic eff'ects, we have also implemented aperture photome- 
try (hereafter AP). For the AP, given an object of certain angular 
size R, the method evaluates the mean temperature in a circle 
of radius r = R and subtracts from it the average found in a 
surrounding ring of inner and outer radii r = R and r = fR, 
respectively, with / > 1 (see e.g., Hernandez-Monteagudo & 
Rubino-Martm 2004). By removing the mean temperature in the 
outer region, the method corrects for large-scale fluctuations in 
the background. Once the temperature estimates are derived for 
each frequency map, they are combined with inverse- variance 
weighting to derive an SZ signal estimate by using the known 
frequency dependence of the (non-relativistic) thermal SZ ef- 
fect. Our choice for the two parameters of the AP method is 
(RJ) = (FWHM, V2). Note that the FWHM varies from one 
frequency band to another. We also note that the flux estimates 
within the aperture have to be corrected separately at each fre- 
quency by an appropriate factor in order to obtain the total 
flux of the source. For example, if the objects are unresolved 
and we assume Gaussian beams, then this correction factor is 
(1 -h exp(-8 In 2) - exp(-4 In 2))"^ for the above choice of R and 
/. For extended objects (e.g., those objects with R500 larger than 
the beam size, and which are modeled here using the "universal 
pressure profile"), the conversion factor can be evaluated numer- 
ically. 

Using one of these methods (MMF or AP), we obtain a 
measure of the intrinsic SZ signal strength ?5oo(0 and the as- 
sociated measurement uncertainty d-e^QQ(i) for the halo of each 
galaxy /. The majority of these individual SZ measurements 
have low signal-to-noise ratio. Following the approach in Planck 
Coflaboration X (2011) and Planck Coflaboration XII (2011), 
we bin them by stellar mass, calculating the bin-average signal 
(Ysooh = [Z Y5oo(i)/d-lJi), with uncertainty 

(T^^ = Zjl^ ^/^^5oo^^^' ^^^^^ is the number of galaxies in bin 
b. 

4. Results 

Our main observational result is given in Fig. 4 and Table 1, 
showing the mean SZ signal measured using the three-band 
MMF for locally brightest galaxies binned according to stel- 
lar mass. In the plot, the thick error bars show the uncertainty 
propagated from the individual measurement errors as described 
above, while the thin bars with large terminators give the vari- 
ance of the weighted bin-average signal found by a bootstrap 
resampling. For the latter, we constructed 1,000 bootstrap real- 
izations of the original LBG catalogue and performed the full 
analysis on each. 

The inset uses a linear scale to better display the significance 
of our detections. We have a clear signal down to the bin at 
11.2 < logio(M,/Mo) < 11.3, centred at M, = 1.8 x lO^^Mo. 
The next three bins provide evidence that the signal continues to 
lower mass with "detections" significant at the 1.6 cr, 1.6 cr and 
2.6 cr levels, from high to low mass, respectively. The last bin 
is centered at = 9 x 10^^ M©, corresponding to a mean halo 
mass of M200 ~ 1.4 X 10^^ M©. These last three bins, however, 
are more seriously aff'ected by dust contamination, as discussed 
below, and for this reason may be more uncertain than these sta- 
tistical measures suggest. 

5. Systematic errors 

In this section, we present a number of tests of the robustness 
of our principal result against systematic error. In the foUow- 
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Fig. 4. Mean SZ signal vs. stellar mass for locally brightest 
galaxies. Thick error bars trace the uncertainty on the bin av- 
erage due purely to measurement error, while thin bars with 
large terminators show the variance calculated by bootstrap re- 
sampling and so also include the intrinsic scatter in the sig- 
nal. The inset provides a view on a linear scale to better eval- 
uate the significance of the detections. We observe a clear re- 
lation between the mean SZ signal and stellar mass down to 
logio(M=,/Mo) = 11.25 (the detection in this bin is at 3.5 cr), 
with a suggestion of signal to lower mass: the next three bins 
show signal at 1.6 cr, 1.6 cr and 2.6 cr, respectively. 

Table 1. Planck SZ signal measurements ?5oo binned by stellar 
mass (adopting a WMAP7 cosmology). These data are displayed 
in Fig. 4. 



^501 



[10- 



Errors [10 ^ arcm^] 
Statistical Bootstrap 



10.05 
10.15 
10.25 
10.35 
10.45 
10.55 
10.65 
10.75 
10.85 
10.95 
11.05 
11.15 
11.25 
11.35 
11.45 
11.55 
11.65 
11.75 
11.85 
11.95 



0.47 


+0.45 


+0.44 


0.79 


+0.41 


+0.39 


0.44 


+0.39 


+0.37 
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+0.37 


0.05 


+0.37 


+0.34 


0.65 


+0.38 


+0.37 


0.80 


+0.39 


+0.40 


0.25 


+0.43 


+0.43 


-0.05 


+0.50 


+0.75 


1.54 


+0.60 


+0.58 


1.27 


+0.78 


+0.78 


1.7 


+ 1.0 


+ 1.1 


5.2 


+ 1.5 


+ 1.8 


11.2 


+2.3 


+2.4 


29.0 


+3.6 


+3.8 


60.7 


+6.2 


+6.8 


123 


+ 11 


+ 16 


266 


+23 


+36 


445 


+53 


+84 


721 


+ 103 


+210 



ing, unless otherwise stated, all results use data at 100, 143, and 
217 GHz only. 

5.1. Stacking real-space reconstructed SZ maps 

According to Fig. 4, the lowest bin at which we have a 
> 3cr detection is the one at logiQ(M^/MQ) = 11.25. As 



a consistency check, and also as an illustration of the fre- 
quency dependence of the detected SZ signal. Fig. 5 shows 
stacked images of central galaxies in six diff'erent mass bins 
of width AlogiQM* = 0.2 centred at logiQ(MJM0) = 
11.05, 11.15, 11.25, 11.35, 11.45, and 11.55. 

The stacked maps are obtained, using equal weights, from a 
(full-sky) SZ map constructed from the Planck 100, 143, 217, 
and 353 GHz maps using a modified internal linear combination 
algorithm (MILCA, Hurler et al. 2010) that has been used for 
other Planck Intermediate Papers (e.g., Planck Collaboration V 
2012; Planck Collaboration X 2012). The well-known internal 
linear combination approach (e.g., Eriksen et al. 2004) searches 
for the linear combination of the input maps that minimises the 
variance of the final reconstructed map while imposing spectral 
constraints. This preserves the thermal SZ signal and removes 
the CMB contamination (using the known spectral signatures of 
the two components) in the final SZ map. The resulting map used 
for this analysis has an angular resolution (FWHM) of 10'. We 
have checked that almost identical maps are obtained with other 
methods. 

The SZ signal is clearly visible in all panels with 
logio(M=,/Mo) > 11.25. The stacked maps show no sign of a 
gradient in the residual signal in the vertical direction, showing 
that the MILCA method is very eff'ective in removing Galactic 
emission. Below the mass limit of 11.25, there is also some ev- 
idence of SZ signal, although here the contrast relative to the 
noise is lower. Finally, we note that the signal in the lower stellar 
mass panels is extended. This is mainly due to the larger satellite 
fraction at these masses that results in a significant contribution 
to the stack from relatively massive halos with centres signifi- 
cantly off'set from the locally brightest galaxy (see Figs. 2, 3 and 
Appendix C). We also discuss in Appendix C the impact of dust 
contamination on these maps. 

5.2. Null tests 

Null tests used to check for systematic errors are shown in Fig. 6. 
Taking the set of MMF filters adapted to each target galaxy, we 
shift their positions on the sky, either with a random displace- 
ment (i.e., by generating a random distribution of new positions 
isotropically distributed outside our Galactic mask) or by shift- 
ing all coordinates one degree in declination, and rerunning our 
analysis. In both cases the result should be zero. The shifted fil- 
ter sets indeed have bin-average SZ signals consistent with zero 
over the entire mass range. 



5.3. Size effects 

Although most of the halos traced by our locally brightest 
galaxies are, according to their inferred R500 values, at most 
marginally resolved by the Planck beams, size eff'ects are not 
negligible, and the full pressure profile has to be used for the flux 
determination. If instead the objects are (incorrectly) assumed to 
be point-like, we find that the flux is underestimated by roughly 
20-30%, although the slope of the Y^qq-M^ scaling relation is 
practically unafl'ected. 



5.4. Photometry comparison 

Figure 7 compares the SZ signal extracted for our LEG sam- 
ple by the two photometry methods described above, namely 
MMF and AP. Here, we compare the total SZ flux from MMF 
(computed as Ysr^^^) with the total flux recovered from the 
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Fig. 5. Equal- weighted stacks of reconstructed SZ maps (i.e., Comptonization parameter maps) for objects in six mass bins centred, 
from left to right and top to bottom, at logio(M,/Mo) = [11.05, 11.15, 11.25, 11.35, 11.45, 11.55]. In all cases, the bin size is taken 
to be 0.2, so the galaxies in two consecutive panels partially overlap. Maps are 2° on a side, with Galactic north at the top. The SZ 
signal traced by the central galaxies is clearly detected in all bins above logio(M=,/Mo) = 11.25. In all panels, the circles indicate 
the FWHM of the data, which corresponds to 10'. 



AP method, after applying the correction factors described in 
Sect. 3. For simplicity, we assume point-like objects for the flux 
extraction in this analysis. This is why the MMF data points 
difl'er from the corresponding points in Fig. 4. For the AP, we 
also compare the nominal four-band analysis with a three-band 
case to illustrate the impact of residual foregrounds on our flux 
estimates. When the 353 GHz channel is included, the AP flux 
estimates at low stellar-mass are biased towards high SZ val- 
ues. This indicates contaminating high-frequency emission as- 
sociated with the sources, presumably dust in the LBGs or their 
satellites. We discuss this issue further in Sect. 5.5. 

The main conclusion is that the two methods, despite their 
difl'erent data processing approaches, produce fully consistent 
results for logiQ(M=,/Mo) ^ 1 1 .25, while the results start to show 
a dependence on the method for stellar masses below that limit. 

5.5. Dust contamination 

The analysis of the last section suggests that our SZ signal es- 
timates may be contaminated by residual dust emission that in- 



creases with frequency and could bias our primary results. To 
evaluate the potential efl'ects, we have performed measurements 
using three difl'erent MMFs, as shown in Fig. 8. The green trian- 
gles and red diamonds represent the results of using all six HFI 
channels or only the lowest three (100, 143, 217 GHz), respec- 
tively. In both cases there is no explicit allowance for a possible 
dust contribution. The blue crosses show results for a modified 
six-band MMF that includes amplitude fits not only to the SZ 
spectrum, but also to a fiducial thermal dust spectrum. 

The three sets of measurements fully agree for the stel- 
lar masses for which we unambiguously detect the SZ signal, 
logiQ M^/Mq > 11.25. This indicates that dust emission does 
not significantly aff'ect our results for these stellar mass bins. 
At lower mass the three-band results and the dust-corrected six- 
band results remain consistent, but the six-band results without 
explicit dust correction are systematically diff'erent. Dust emis- 
sion is clearly sufficient to contaminate our six-band filter esti- 
mates of SZ signal if uncorrected, but it does not appear to be a 
major problem when only the lower three frequency bands are 
used. The residual dust contribution estimated from the scatter 
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Fig. 6. Null tests performed on the locally brightest galaxy sam- 
ple. Red points correspond to placing the filter one degree in dec- 
lination away from the position of each LBG, while green points 
correspond to random high latitude filter positions. Both sets are 
consistent with zero. The black points show our measurements 
with filters centred on the LBG sample, demonstrating highly 
significant detections. 
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Fig. 7. Comparison of the SZ measurements on the full locally 
brightest galaxy sample for two diff'erent photometry methods: 
the matched multi- filter (MMF) and the aperture photometry 
(AP) approach. For this figure and only for this figure, we as- 
sume point-like objects for both methods and plot the derived 
total SZ flux (or the flux within 57?5oo for MMF). The signal 
detected by the two methods is consistent at all stellar masses 
when only three frequencies are used, but when four frequen- 
cies are used, the AP results are contaminated by high frequency 
emission at stellar masses below -2x10^^ M^. 



and ofl'set of the red and blue points for log^o M^/Mq < 1 1.0 is 
below ~ 10"^ arcmin^ and so lies comfortably below our mea- 
sured signal. 

There is a clear indication of signal in the three bins just be- 
low logiQ M^/Mq = 1 1.25 both for the three-band MMF and for 
the dust-corrected six-band MMF. However, the dust-corrected 
results appear systematically lower than the (uncorrected) three- 
band results by an amount similar to that seen at lower masses 
where the SZ signal is undetected. Further, the six-band MMF 
measurements without dust correction (the green triangles) dif- 
fer substantially for these (and all lower) bins. This suggests that 
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Fig. 8. Impact of dust contamination on our SZ measurements. 
Three cases are shown: a 6-band MMF (all Planck HFI frequen- 
cies) with no explicit allowance for a dust contribution (green 
triangles), a 3 -band MMF also with no explicit dust modelling 
(red diamonds); and a modified 6-band MMF that includes an 
amplitude fit to a fiducial dust spectrum (blue crosses). The error 
bars include measurement uncertainties only. For stellar masses 
where we clearly detect the signal (i.e., at log^QM^^/M© > 
1 1.25), the three measurements agree, indicating that dust emis- 
sion does not significantly aff'ect those measurements. At lower 
masses the 3 -band results are consistent with the 6-band results 
when dust is explicitly included in the modelling, but not other- 
wise. 



dust emission aff'ects these stellar mass bins noticeably even for 
the three-band MMF, so the corresponding points in Fig. 4 may 
be more uncertain than indicated by their statistical error bars. 
Although formally the dust-corrected six-band MMF would ap- 
pear to give our most accurate estimates of stacked SZ signal, 
we are uncertain whether the fiducial dust spectrum it assumes 
is appropriate for these specific sources. Therefore we conser- 
vatively quote results based on the three-band MMF, using the 
dust-corrected six-band results to give an estimate of remaining 
dust-related systematics. 

Finally, we note that residual dust contamination bi- 
ases the (uncorrected) six-band MMF signal estimates for 
logiQ M^/Mq < 11 (Fig. 8) in the opposite direction to the AP 
signal estimates (see Fig. 7). The agreement of the two meth- 
ods for logiQ M^/Mq > 11.25 is thus a further indication of the 
robustness of our primary results. 

5.6. Stability of the signal in different sky surveys 

We have also checked that the SZ signal is stable against split- 
ting the Planck data into complementary subsets. For instance, 
the signal obtained from the maps of the first 6 months of ob- 
serving time is fully consistent with that obtained from maps of 
the second 6 months and the last 3.5 months (of course the latter 
has larger error bars due to its smaller sky coverage). 



6. The F500-M500 relation 

We now turn to the interpretation of our measurements in terms 
of the SZ signal-halo mass scaling relation: Y500-M500. Our con- 
clusions are summarised in Fig. 9. 

From our simulation of the locally brightest galaxy cata- 
logue, we expect a large range of halo masses within a given 
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Fig. 9. Left: Comparison of the measured mean SZ signal as a function of LBG stellar mass (red points) to simulated observations 
(blue points). The simulations assign to each observed LBG the halo mass and positional offset of a randomly chosen simulated 
LBG of the same stellar mass (compare Fig. 3). Our best fit 7500-^500 scaling relation is then used, together with the universal 
pressure profile, to inject a simulated signal into the Planck maps (see text). An "observed" signal is obtained by applying the MMF 
exactly as for the real data. The inset gives the ratio of the bin-averaged injected and actual signals. Right: Mean SZ signal as a 
function of eff'ective halo mass. The bin-averaged SZ signal measurements of the left panel have been translated to this plane using 
the simulations as described in the text (the red points). The dot-dashed line is our best fit relation between halo mass and SZ signal, 
i.e., the one leading to the simulated measurements in the left panel. The green points give the mean SZ signal of MCXC clusters 
binned by a halo mass estimated from their X-ray luminosity using the REXCESS relation without correction for Malmquist bias 
(line 3 in Table 2 of Planck Collaboration X 201 1). The dashed blue line shows the self-similar model calibrated on the REXCESS 
sample as given by Arnaud et al. (2010). The inset gives the ratio of all measurements to this model's predictions. As in previous 
figures, the thick error bars account only for measurement uncertainties, while thin bars with large terminators result from a bootstrap 
analysis and so include intrinsic scatter eff'ects. 



bin of stellar mass and, in addition, a fraction of galaxies that 
are, in fact, satellites, with significant positional off'sets relative 
to their host halo (see Fig. 3). These eff'ects impact our measure- 
ments of the SZ signal- stellar mass relation in two ways. First, 
the MMF is not perfectly matched to each individual object be- 
cause we fix the filter scale to the median halo size. This causes 
an aperture-induced bias in the flux measurement. Second, our 
filter is miscentred for those systems where the LBG is, in fact, 
a satellite. These galaxies are often associated with substantially 
more massive dark halos than typical LBGs of the same stel- 
lar mass, leading to an increase in the mean signal in the bin, 
mitigated by the substantial angular ofl'sets of most such satel- 
lites from the true centres of their rich clusters. This increases 
the apparent extent of the signal in stacked maps like Fig. 5, but 
decreases the contribution to the signal through a matched filter 
centred on the galaxy (see Appendix C). 

Using our simulation of the LBG catalogue, we can account 
fully for these eff'ects and extract the underlying Fsoo-^soo rela- 
tion in an unbiased way. Within each stellar mass bin, we iden- 
tify each observed LBG with a randomly chosen simulated LBG 
of the same stellar mass, assigning it the halo mass and positional 
off'set from halo centre of its partner, but retaining its observed 
redshift. We give each halo a SZ signal distributed according to 
the "universal pressure profile" and normalized using a specific 
model 7500-^500 scaling relation. Each synthesised object is then 
observed with the three-band MMF centred on the galaxy's posi- 
tion, and the measurements are binned and weighted in the same 
way as the real data to obtain (F)s. 

This procedure enables us to translate a model 7500-M500 re- 
lation to our observational plane, Fsoo-M.,, and thus to fit for the 



underlying scaling relation with halo mass M500. We model this 
relation as 



?500 = 



M500 



3xl0i4M(: 



(1) 



fixing the mass exponent to its self-similar value, = 5/3, 
and fitting for the normalization 1m • Restricting the fit to 
XogiQ^M^IMo) ^ 11.5, for direct comparison to X-ray samples 
in the discussion below, we find 



Ym = (0.73 ± 0.07) X 10"^ arcminl 



(2) 



In the left-hand panel of Fig. 9, the red points reproduce the mea- 
surements given in Fig. 4, while the blue points show the simu- 
lated observations for this best-fit ^500-^500 scaling relation. 

The best-fit is, however, formally unacceptable, with a re- 
duced of 3, which we can more readily appreciate from the 
inset showing the ratio of the actual observations to the simu- 
lated bin averages on a linear scale. The data prefer a shallower 
slope than the self-similar o^m = 5/3 over the mass range of the 
fit. Moreover, we see that a power law cannot fit the data over 
the full mass range probed by our measurements. To ease com- 
parison with the X-ray sample, we will nevertheless adopt this 
fit below. 

The right-hand panel of Fig. 9 considers the SZ signal-halo 
mass plane. The blue dot-dashed line simply traces our best- fit 
Y500-M500 relation. The blue dashed line is the self- similar re- 
lation derived from X-ray cluster studies (Arnaud et al. 2010), 
while the green points present binned SZ measurements for the 
approximately 1,600 clusters in the Meta-Catalogue of X-ray de- 
tected Clusters (MCXC) (Piff'aretti et al. 2010). The latter mea- 
surements are as reported in Planck Collaboration X (201 1), with 
one minor change: in Planck Collaboration X (2011) we used 
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an empirical slope for the Y500-M500 relation taken from X-ray 
studies; for the points in Fig. 9, we repeated the same analy- 
sis fixing the slope instead to its self-similar value, as was done 
for the LBG sample. This change moves the green points only 
very slightly relative to those shown in Planck Collaboration 
X (2011). For the mass estimates of the MCXC objects, we 
applied the X-ray luminosity-mass relation from Pratt et al. 
(2009), corresponding to the case of line 3 of Table 2 in Planck 
Collaboration X (201 1). The mass is calculated for each MCXC 
cluster and then binned. We plot the point at the median value of 
the mass in each bin. 

To transcribe our central galaxy catalogue measurements 
onto this figure, we must first find the eff'ective halo mass corre- 
sponding to each stellar mass bin. This eff'ective mass is a com- 
plicated average over the halo masses within the bin, weight- 
ing by the fraction of SZ signal actually observed, i.e., after ac- 
counting for aperture and miscentering eff'ects. The bin-averaged 
mean SZ signals we estimate for our mock LBG catalogue in- 
clude all these eff'ects, and so can be used to calculate an eff'ective 
mass as M^^ = 3 x 10^^ ({Y}JYm)^^'''^, where {¥% is calcu- 
lated for each bin as described above, and Ym and o^m = 5/3 are 
the parameters used for Eq. 1 in the simulation. (Note that the 
result is independent of the normalisation Ym-) We do this for a 
suite of simulated catalogues and take the ensemble average ef- 
fective mass for each bin, plotting the results as the red points in 
the right-hand panel of the figure. 

These LBG results extend the SZ-halo mass scaling relation 
down in mass by at least a factor of 3, to M500 = 2 x 10^^ M© 
(the stellar mass bin at log^Q M^/Mq = 11.25). This is the low- 
est halo mass for which the mean SZ signal has been measured. 
As previously discussed, there is a clear indication that the re- 
lation continues to even lower mass, with marginally significant 
detections in the next three stellar mass bins. The lowest stellar 
mass bin with an apparent SZ detection (at 2.6cr) corresponds 
to eff'ective halo mass log^o Msoo/M© = 12.6. Our power-law fit 
adequately describes the data points over more than two orders 
of magnitude in halo mass down to this remarkably low value 
with no hint of a significant deviation. 

The inset in the right panel of Fig. 9 shows the ratio of our 
measured mean SZ signal to that predicted by the self- similar 
scaling relation deduced from X-ray observations of clusters 
(the dashed blue line (Arnaud et al. 2010)). Direct measure- 
ments obtained by binning the MCXC clusters (the green points) 
agree with this relation. This was the principal result of Planck 
Collaboration X (2011). The SZ measurements for our LBGs 
fall below the relation, however. The horizontal dot-dashed line 
gives the ratio our LBG fit to the X-ray model (this is the off'set 
between the two blue lines in the main figure). Recall that the fit 
to the LBG catalogue was restricted to masses overlapping the 
X-ray sample, log^g M5OO/M0 > 13.8. Over this range, the mean 
SZ signals associated with LBG halos are about 20 % lower than 
found for X-ray clusters with the same halo mass, a diff'erence 
that is significant at the 2.6 cr level. 

A number of eff'ects could contribute to an off'set of this 
size. The masses plotted for the MCXC were calculated using a 
luminosity-mass relation derived from the REXCESS sample as- 
suming that halo mass scales self- similarly with the mass-proxy 
Yx and without correction for Malmquist bias (Pratt et al. 2009). 
Using the Malmquist-corrected relation would remove much of 
the off'set and bring the two F500-M500 scaling relations into ac- 
ceptable agreement. In this sense, the off'set is consistent with 
the estimated eff'ects of Malmquist bias on the X-ray sample. 
However, such biases depend on the detailed selection proce- 
dure of the stacked and calibrating cluster samples, on the way 



in which the calibration relation is derived, and on the (corre- 
lated) intrinsic scatter of clusters around the Lx-Msqq and 7500- 
M500 relations. Thus they can only be corrected through detailed 
modelling both of the cluster population itself and of the defi- 
nition and analysis of the specific cluster surveys involved (e.g., 
Angulo et al. 2012). Furthermore, halo masses are estimated in 
very diff'erent ways in our two samples — from X-ray luminosi- 
ties calibrated against individual hydrostatic mass measurements 
for the MCXC, and through an abundance matching argument 
based on the WMAP7 cosmology for the LBG catalogue. Any 
off'set between these two halo mass scales will result in off'sets 
in Fig. 9. For example, a number of recent papers have argued 
that failure of some of the assumptions underlying the standard 
methods for estimating cluster masses from X-ray data (e.g., de- 
tailed hydrostatic equilibrium or the unimportance of turbulent 
and nonthermal pressure) could produce a systematic bias in the 
X-ray cluster mass scale (Planck Collaboration XII 2011; Rozo 
et al. 2012; Sehgal et al. 2012). Finally, as for the LBG sam- 
ple, each luminosity bin of the MCXC contains a distribution of 
halo properties that are averaged in complicated fashion by our 
stacked SZ measurement. Understanding the relative importance 
of these various eff'ects at a precision better than 20% would 
again require detailed modeling of the heterogeneous MCXC 
catalogue. 

7. Conclusions 

Using Planck data, we have measured the scaling relation be- 
tween Sunyaev-Zeldovich signal and stellar mass for locally 
brightest galaxies (Ysqq-M^). This is the first time such a rela- 
tion has been determined, and it demonstrates the presence of 
hot, diff'use gas in halos hosting central galaxies of stellar mass 
as low as M:, = 2 X 10^^ M©, with a strong indication of signal at 
even lower masses. We have constructed a large mock catalogue 
of locally brightest galaxies from the Millennium Simulation 
and used it to model the Planck observational process in detail 
in order to extract from our measurements the underlying SZ 
signal-halo mass relation (Ysoq-Msoq). This new relation spans 
a large range in halo mass, reaching from rich clusters down to 
M500 = 2.0 X 10^^ Mq, with a clear indication of continuation to 
^500 - 4 x 10^^ Mq. This is the lowest mass scale to which an 
SZ scaling relation has so far been measured. The fact that the 
signal is close to the self-similar prediction implies that Planck- 
detected hot gas represents roughly the mean cosmic fraction 
of the mass even in such low-mass systems. Consistency with 
their low observed X-ray luminosities then requires the gas to be 
less concentrated than in more massive systems. Integration of 
the halo mass function down to M500 = 4 x lO^^M© shows that 
Planck has now seen about a quarter of all cosmic baryons in the 
form of hot gas, about four times as many as are inferred from 
X-ray data in clusters with M500 > lO^'^M©. 

At the high mass end, the scaling relation we derive from 
our LBG data shows reasonable agreement with X-ray cluster 
results. The 20% lower normalisation that we find (significant 
at the 2.6(T level) can be explained in principle by a number of 
possible eff'ects related to the diff'ering selection and mass esti- 
mation methods of the two samples. Agreement at this level of 
precision is remarkable, and understanding the remaining diff'er- 
ence would require detailed modeling of the selection and cali- 
bration of the X-ray samples. The fact that plausible Malmquist 
corrections can eliminate most of the diff'erence shows that clus- 
ter studies are now reaching the ~ 10% precision level. 

We find that the 7500-^500 scaling law is described by a 
power law with no evidence of deviation over more than two 
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orders of magnitude in halo mass. The gas properties of dark 
matter halos appear remarkably regular over a mass range where 
cooling and feedback processes are expected to vary strongly. 
In particular, we find no change in behaviour in the low-mass 
systems for which substantial feedback eff'ects are invoked in 
current galaxy formation models (e.g., from AGN). Statistical 
studies of large galaxy and cluster samples, such as those pre- 
sented here, can clearly shed new light on the thermal cycle at 
the heart of the galaxy formation process. 
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Appendix A: Robustness of our results to 
variations in isolation criteria 

As explained in Sect. 2.2, our locally brightest galaxy catalogue 
was built starting from a parent population with r < 17.7 taken 
from the spectroscopic NYU-VAGC and eliminating any can- 
didate with a companion of equal or brighter r magnitude vio- 
lating certain isolation criteria. In particular, we defined locally 
brightest galaxies to be the set of all objects with z > 0.03 that 
are brighter than all other sample galaxies projected within a ra- 
dius of Riso = 1.0 Mpc, and diff'ering in redshift by less than 
1,000 km s-i. Hereafter, we refer to these criteria as the "1 Mpc 
case". 

To test the robustness of our results against changes in these 
isolation criteria, we compared them to a case with stricter iso- 
lation criteria, Ri^o = 2.0 Mpc and 2,000 km s"^ in redshift, here- 
after the "2 Mpc case". 

Applying the isolation criteria to the parent spectroscopic 
catalogue as before, but with these new values, we end up with 
a first sample of 206,562 locally brightest galaxes. Again, in a 
second step we use SDSS photometry to further eliminate ob- 
jects with companions that might violate the isolation criteria. 
After removing any candidate with a (photometric) companion 
of equal or brighter r magnitude and projected within 2.0 Mpc, 
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we end up with a cleaned sample of 110,437 locally brightest 
galaxies. In particular, this sample contains 58,105 galaxies sat- 
isfying the bound log^o^^/M© > 11.0, which is the regime 
where we find significant SZ signal. Thus, 23,287 galaxies in 
this mass range are eliminated from the sample studied in the 
main body of this paper by the stricter isolation criteria. 

To evaluate the reliability of the new Ri^o = 2.0 Mpc sample, 
we follow the same procedure as before and construct a mock 
sample of locally brightest galaxies from the Guo et al. (2011) 
simulation. As expected, the new Ri^o = 2.0 Mpc mock sample 
has a higher reliability than the Ri^o = 1.0 Mpc case. The fraction 
of locally brightest galaxies that are centrals now has a minimum 
of just over 87 % at stellar masses somewhat above 10^^ M©. The 
improvement is less than might have been anticipated because, 
as noted in Sect. 2.2. 1 , the majority of the satellite galaxies in our 
simulated 1 Mpc sample were included because they are brighter 
than the central galaxies of their own halos, rather than because 
the isolation criteria failed to eliminate them. 

Finally, Fig. A.l compares the SZ signal-halo mass relations 
(Y500-M500) derived for the two cases (1 Mpc and 2 Mpc). Halo 
masses for the 2 Mpc case are computed as explained in Sect. 6 
(see Table B.l for the numerical values). The main conclusion is 
that the SZ signal-halo mass scaling relation is not sensitive to 
the isolation criteria. 



Appendix B: Predicted properties of the stellar 
mass-halo mass relation 

Using our mock catalogues based on the semi-analytic galaxy 
formation simulation of Guo et al. (2011), we provide here ad- 
ditional information on the predicted properties of the stellar 
mass-halo mass relation. Figure B.l shows the distribution of 
halo mass (Mh) predicted for nine of the stellar mass bins con- 
sidered in this paper, and for two sets of isolation criteria: the 
1 Mpc and 2 Mpc cases (see Sect. 2.2 and Appendix A). Vertical 
lines correspond to the mean (red), median (green), and the "ef- 
fective" (blue) values of halo mass in each bin. The correspond- 
ing numbers are listed in Table B.l, which also gives the RMS 
of the posterior Mh distribution. The eff'ective halo masses are 
computed as described in Sect. 6. 



Appendix C: Impact of miscentering and scatter on 
the binned SZ signal and stacked SZ maps 

As discussed in Sec. 6, we used the semi-analytic galaxy forma- 
tion simulation of Guo et al. (201 1) to account for the eff'ects of 
miscentering and scatter in halo mass at fixed stellar mass when 
interpreting our measurement (see Figs. 4 and 9). Figure C.l 
isolates the impact of each eff'ect on the binned SZ measure- 
ments, using the procedure outlined in that section. The green 
points represent the ideal case with no miscentering and SZ filter 
perfectly matched to the size of each individual object. The red 
crosses add miscentering off'sets taken from the off'set distribu- 
tion in the simulations for each stellar mass bin. The drop in SZ 
amplitude is expected because we now miss SZ signal from the 
miscentered objects. Additionally fixing the filter size accord- 
ing to the median halo mass in each stellar mass bin, as done 
throughout this paper, we recover our previous results, shown as 
the blue triangles here and as the red diamonds in the left-hand 
panel of Fig. 9. 

Using the same simulations, we can also estimate the im- 
pact of miscentering on the stacked SZ maps of locally brightest 



galaxies (see Fig. 5). Here, we use the full simulation to com- 
pute Tp, the projected distance of each locally brightest galaxy 
from the gravitational potential minimum of its halo. Average 
and RMS values for rp for all the stellar mass bins considered 
in this paper and for the 1 Mpc and 2 Mpc samples are given in 
Table C.l. Histograms of these rp values are shown in Fig. C.2. 
Note that the median value of rp, which is not listed in the table, 
is zero for all bins. 

These values can be used to predict the impact of miscenter- 
ing of the locally brightest galaxy with respect to its halo (and 
thus, with respect to the centre of the associated SZ emission). 
Figure C.3 illustrates the broadening of the SZ stacked profile 
caused by this eff'ect. For this computation, we assume point-like 
objects and a Gaussian beam profile of 10' for easier compari- 
son with Fig. 5. For each stellar mass bin, the Mh value from the 
simulation is used to predict the total SZ flux using Eqs. 1 and 2, 
and the rp value is used to ofl'set the position of the SZ signal. In 
order to convert rp values (in physical units) into angular ofl'sets, 
a redshift for each simulated object is drawn from the observed 
distribution for locally brightest galaxies of similar stellar mass. 
Miscentering broadens the stacked SZ profile, yielding typical 
FWHM of ~ 20' for log^^MJUe < 11.25, and also modifies 
the shape of the profile, by increasing the amount of SZ flux in 
the tails of the distribution. These values are slightly smaller (but 
comparable) to the observed widths of the SZ emission in Fig. 5. 

Finally, Fig. C.4 shows equal-weighted stacks of SZ maps 
centred on the real central galaxy sample, similar to those of 
Fig. 5, but now using all six HFI frequency channels in the 
MILCA algorithm, rather than just the lowest four. For all six 
stellar mass bins the noise in these new maps, as measured by 
the RMS fluctuation about the mean in pixels more than 20' from 
map centre, is lower than in the maps of Fig. 5. This shows that 
the addition of high frequency information has improved the ac- 
curacy with which non-SZ signals, primarily dust emission, are 
removed. Almost all this improvement comes from the inclu- 
sion of the 545 GHz channel; maps made with and without the 
857 GHz channel are almost identical. As a result of this im- 
provement, the signal-to-noise ratio of the peaks near the map 
centre is higher in all the panels of Fig. C.4 than in the corre- 
sponding panels of Fig. 5. This strengthens our conclusion that 
the apparent SZ signals near the centres of the two lowest stellar 
mass panels are, in fact, real, despite their apparent breadth and 
irregularity. The breadth is likely due to the miscentering efl'ects 
explored above while the irregularity looks consistent with the 
overall noise level of the maps. 
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Fig. A.l. Left: Comparison of the SZ signal-halo mass scaling relation for two different sets of isolation criteria. The triple-dot 
dashed line is our best fit model (see Eq. 1 and 2). Right: Same as above, but now showing the ratio of the previous data points to 
the Arnaud et al. (2010) Y-Msm relation. 
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Fig.B.l. Probability distribution function of halo mass, Mh, for nine of the stellar mass bins considered in this paper. Solid lines 
correspond to the sample isolated according to the 1 Mpc criteria, while dashed lines show the distributions for the 2 Mpc sample. 
Vertical colored lines show three diff'erent characteristic masses (the mean, median, and "eff'ective" halo masses) for the 1 Mpc 
sample (see Table B.l for numerical values). 
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Table B.l. Statistics of halo mass for various stellar mass bins, for the 1 Mpc and 2 Mpc isolation criteria. The first three columns 
for each case (mean, median, and RMS values for the halo mass) are derived from the simulation only, while the eff'ective halo mass 
uses the redshifts and stellar masses of the observed galaxies, as described in Sect. 6. All masses (M) in this table are decimal 
logarithms of the value in units of M©. 
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Fig. C.l. Impact of miscentering and scatter on the binned SZ 
measurements. The green points give results in an ideal situation 
with no miscentering and SZ filter perfectly matched to each in- 
dividual object in a given stellar mass bin. The red crosses add 
the effect of miscentering, with offsets drawn from the distribu- 
tions given by the simulations for each stellar mass bin. The blue 
triangles additionally include the aperture effect caused by fixing 
the filter size according to the median value of the halo mass in 
each bin. 
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Table C.l. Statistics of the distribution of distances of the lo- 
cally brightest galaxies from the gravitational potential minima 
of their parent halos, for the 1 Mpc and 2 Mpc isolation criteria. 
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Fig. C.2. Distribution of offsets of locally brightest galaxies from the gravitational potential minima of their parent halos, both for 
the 1 Mpc (black) and for the 2 Mpc (green) isolation criteria. Table C.l gives mean and RMS values for these distributions. 
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Fig. C.3. Impact of miscentering on stacked SZ maps. See the text for details of the simulation shown here. For an original resolution 
of FWHM= 10', miscentering broadens the stacked profiles to a FWHM- 20' for log^o MJMq < 11.25. 
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Fig. C.4. Similar to Fig. 5, but using a reconstructed SZ map that now uses all six HFI frequency channels. The noise in all maps is 
reduced by the inclusion of the two highest frequencies. Stacked images in the stellar-mass bins above logiQCM^^/M©) = 11.25 are 
not significantly aff'ected, but for the low stellar-mass panels, the extended signal near map centre is larger and has higher signal to 
noise than in Fig. 5, suggesting that it may be real SZ signal broadened by miscentering eff'ects. 
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